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Abstract. We present a detailed study of the finite one-dimensional quantum Ising 
chain in a transverse field in the presence of boundary magnetic fields coupled with 
the order-parameter spin operator. We consider two magnetic fields located at the 
boundaries of the chain that have the same strength and that are aligned in the same 
or in the opposite direction. We derive analytic expressions for the gap in all phases for 
large values of the chain length L, as a function of the boundary field strength. We also 
investigate the behavior of the chain in the quantum ferromagnetic phase for oppositely 
aligned fields, focusing on the magnet-to-kink transition that occurs at a finite value of 
the magnetic field strength. At this transition we compute analytically the finite-size 
crossover functions for the gap, the magnetization profile, the two-point correlation 
function, and the density of fermionic modes. As the magnet-to-kink transition is 
equivalent to the wetting transition in two-dimensional classical Ising models, our 
results provide new analytic predictions for the finite-size behavior of Ising systems in 
a strip geometry at this transition. 
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The quantum Ising chain is a useful theoretical laboratory in which fundamental issues 
concerning quantum many-body systems can be thoroughly investigated, exploiting the 
exact knowledge of several features of its phase diagram and quantum correlations. 
Many results for its low-energy properties have been derived in the quantum ordered 
and disordered phases, and, in particular, at the quantum critical point separating the 
two phases, both in the thermodynamic limit and in the finite-size scaling (FSS) limit 
with several types of boundary conditions, see, e.g., Refs. [H 2j [3], SJ [5j |6, .7] ® |9j TUI ITT- 
121 EH HU EH EH El HH1 HHJ ED] EH [221 [23] and references therein. 

In this paper we extend these analytic results to the case in which boundary fields 
are present. We present a detailed study of the quantum Ising chain in a transverse 
magnetic field |5] in the presence of magnetic fields coupled with the order-parameter 
spin operator, located at the boundaries of the chain. We assume the two magnetic fields 
to have the same strength and consider two cases: (i) the two fields are parallel; (ii) 
the two fields are oppositely aligned. As expected, in the quantum paramagnetic phase 
the boundary fields do not change the large-size behavior of low-energies quantities. At 
the critical point, bulk behavior is independent of boundary conditions. However, the 
magnetic fields induce a surface phase transition with a corresponding scaling behavior. 
The quantum ferromagnetic phase is more interesting. If the boundary fields have 
opposite direction, one observes two different bulk phases. For small magnetic fields, 
the ground state is ferromagnetic as it occurs in the absence of boundary interactions. 
On the other hand, kink states [1J are the relevant low-energy excitations for large 
boundary fields. The two different phases are separated by a continuous transition that 
is only characterized by the nature of the two coexisting phases and, indeed, the same 
transition occurs in Ising rings in the presence of a localized link defect [23], 

Because of the quantum-to-classical mapping, our results can also be applied 
to the two-dimensional Ising model, and, more generally, to any model in the Ising 
universality class, in a strip geometry. In the two-dimensional case the quantum 
transition between the ferromagnetic and the kink phase corresponds to the wetting 
transition, which, in the context of a strip geometry, is sometimes identified as an 
interface localization/delocalization transition [2T SSI ES, 27]. Our results, therefore, 
provide new analytic expressions for two-dimensional Ising systems at the wetting 
transition in the presence of boundary magnetic fields. 

The paper is organized as follows. In Sec. [2] we introduce the one-dimensional 
quantum Ising chain with boundary fields. In Sec. [3] we compute the low- 
energy spectrum by exploiting the equivalent quadratic fermionic formulation of the 
Hamiltonian Explicit calculations are reported in the following sections. Secs. E 

[5] and [6] report results for the quantum paramagnetic phase, at the quantum critical point 
and for the ordered magnetized phase, respectively. Sec. [7] is devoted to the study of the 
magnet-to-kink transition driven by the boundary fields in the ordered phase. We obtain 
exact results for the gap, the magnetization profile, the two-point spin-spin correlation 
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function, and the entanglement entropy. In Sec. 17.51 these results are compared with the 
existing ones for two-dimensional classical models. Finally, in Sec. [8] we summarize the 
main results and draw our conclusions. A number of appendices report the derivations 
of some of the results. 


2. Model and definitions 


The Hamiltonian of the quantum Ising chain in a transverse magnetic field is given by 




o. 


(3) 


( 1 ) 


i —1 


i =1 


where a O are the Pauli matrices. In the following we assume ferromagnetic nearest- 
neighbor interactions with J — 1, and g > 0. 

The Ising chain undergoes a continuous transition at g = 1 [5], separating a 
quantum ferromagnetic phase (g < 1) from a quantum paramagnetic phase (g > 1). 
In this paper we investigate the effects of boundary magnetic fields aligned along the x 
axis. They give rise to an additional energy term 


H b 


-Cl O’ 


(i) 

i 



( 2 ) 


that is added to Hamiltonian fll]). In the following we shall consider the model in the 
specific case of equal parallel boundary fields (PBF) 


Cl = Ci = C, (3) 

and of equal oppositely-aligned boundary fields (OBF) 

Cl = -Ci = C (4) 


It is not restrictive to assume C > 0 in both cases. 

We will often use a basis in which cC^ is diagonal. States will be labelled as 
|si, S 2 , • ■ • sl), where s* is the eigenvalues of Signs will be fixed so that crf^ has the 
form 



0 1 
1 0 


(5) 


in this basis. In the absence of boundary fields, the Hamiltonian commutes with the 
generator P z = cq- 3 ' ) of the Z 2 transformations |s) —>■ — |s) and with the operator T 

of the reflection transformations defined by T|si, s 2 ,... sfi) = |sl, sl- i, ..., si). These 
operators do not generally commute with H b . However, note that [P Z T, H b \ = 0 for 
OBF. 

In this paper we analyze the low-energy spectrum of the model. In particular, we 
compute the energy differences between the lowest states and the ground state 


A = Ax, 


A n — En Eq, 


( 6 ) 
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(here E n are the energy eigenvalues ordered so that E 0 < Ei < E 2 ..., and A is the 
gap), the local magnetization and the two-point correlation function, 


m{i) = (or, a) ), G(i,j) = (afi'a) 1 ’). 


- /AbRiL 


(7) 


For OBF, because of the symmetry under P 3 T, we have m{i ) = — m(L — i) and, as a 
consequence, the average magnetization m (i) always vanishes. We also define the 
integrated correlation x and the correlation length £ with respect to the center of the 
chain: 

X = E G ( L / 2 ’*), ? = ^ - V2) 2 G(L/2, i). (8) 


3. Jordan-Wigner representation and Hamiltonian diagonalization 


3.1. Fermionic representation 


To determine the spectrum of Hamiltonian ([TJ), we extend the model, considering two 
additional spins located in 0 and L + 1 and the Hamiltonian 

L—1 L 




cr 


(i)Ri) 


a i +1 JO Cq 




- J L aya 


'L l 


(i) 

L +1 




o, 


(3) 


(9) 


i =1 


%— 1 


This is the Ising Hamiltonian with two different couplings on the boundary links and 
zero transverse field on the boundaries. Let us note that the Hamiltonian H e commutes 
with both ctq 1 ' and which can therefore be simultaneously diagonalized. The 

Hilbert space can be divided into four sectors, which we label as (1,1), (—1,1), (1, —1) 
and (—1, —1), where (s 0 ,sl+i) are the eigenvalues of cr^ and crj*^. The restriction of 
H e to each sector gives rise to the Hamiltonian H, defined in Eq. (JTj) , with a boundary 
term of the form (J2j). Hence, the spectrum of H e also provides the spectrum of H + H^. 
Let us also note that H e is Z 2 symmetric. Indeed, if P z = ntJo 1 °i 3) > then i H e, p z] = 0. 
Since P z does not commute with and o^ 1 ^, the spectrum is necessarily degenerate. 
Moreover, since P z anticommutes with these two boundary operators, P z maps sector 
(1,1) to (—1, —1) and (1, —1) to (—1,1), so that the restriction of H e to (1, ±1) allows 
us to compute the full spectrum of H e . 

To compute the spectrum of Hamiltonian ()9|l , we follow Ref. [3]. We first perform 
a Jordan-Wigner transformation, defining fermionic operators c t and c\ 

2—1 

cj = R,afi, Ci = a~R u R t = (—l)* -1 af\ (10) 

3 = 1 

where a ± = (cr^ 1 ) Pa^)/ 2. These relations can be inverted, obtaining 

= Rfic\ + q), af ] = 2 c| Ci - 1. (11) 

Thus, Hamiltonian (JH]) becomes 

L L 

He = ~g ^( 2C ! C * - !) - Ji(c\ C i+ 1 + c\ +1 Ci + cJ C f +1 + Ci+iCi), (12) 
2=1 2=0 
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with J* = J = 1 for i = — 1. In this formalism cr/' 1 = Co + Cq and 

AX\ = P 2 (c^ +1 — Cl + i). To go further, let us rewrite H e as 


L +1 r 


H e = Lg~Y, 


i,j =0 


c l-^u c j d" 2 


4 5 o c ] + 


:CiBj 


-i^ij c j 


(13) 


with A and B symmetric and antisymmetric matrices, respectively. Then, we perform 
a Bogoliubov transformation. We introduce new canonical fermionic variables 

L +1 

flk = + h ki Ci), (14) 

i=0 

where g k % and h kl are fixed by the requirement that H takes the form 


L +1 

H e Egs T ^ ^ ^khkhk) 
k =0 


(15) 


with 0 < To < £\ < ... Following Ref. [3], we define the vectors 

U k = {gw + hko, g k i + Pci, ■ ■ ■)> kfc = ((/fco — Pco, fl'fci — Pci, • • •)• (16) 

The variables rj k satisfy canonical anticommutation relations if the vectors Uk form an 
orthonormal basis, and so does the set V k . The vectors V k satisfy 


(A + B)(A-B)V k = £lV k , 


(17) 


which turns the determination of the energies £ k into an eigenvalue problem. If £ k does 
not vanish, Uk is given by 

U k = X-{A-B)V k . (18) 

If £ k is zero, U k is the null eigenvector of {A — B){A + B). It is also possible to evaluate 
the constant E gs in Eq. (fT5lh which provides the energy of the ground state: 


L -\-1 

Egs = — — £ k - 

k =0 


(19) 


The matrix C = j{A + B)(A — B) is given by (we write the 5x5 matrix C for L = 3, 
the generalization to any L being obvious) 


C 


( Jq gJo 0 0 0 ^ 

g,J 0 1 + g 2 g 0 0 

0 g 1 + g 2 g 0 

0 0 g J 2 L + g 2 0 

\ 0 0 0 0 0/ 


( 20 ) 


This matrix C has a zero eigenvalue, £ 0 = 0, with eigenvector V 0 = (0, ...,0,1). 
Correspondingly, we obtain Uq — (1, 0,..., 0) and 

7o = l(co + 4) + 1(4+1 - Ci+i) = 4 1 ’ + (-l^/W+i 


( 21 ) 
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which depends only on the boundary fermion operators. Note that it can be rewritten as 
a combination of , and P z , hence it does not represent an additional symmetry 

of the Hamiltonian. The presence of a zero eigenvalue was expected, as the spectrum 
is degenerate. If Jo an d Jl do not vanish, since there is only one fermion operator with 
zero energy, the spectrum is doubly degenerate. 

Let us now consider the nonzero eigenvalues. Since the vectors 14 are an 
orthonormal set, orthogonality with Vo implies 14 = (a 0 ,..., a^, 0 ). ft follows 
analogously L4 = (0, hi,, Vl+i). Therefore, we have for k — 1,..., L + 1 

Vk = y ( c o - c o) + _ y^(dL +i + 4 +1 ) + ... ( 22 ) 

where the dots represent a polynomial in Cj and e) with 1 <3 <L. It is then easy to 
verify that 

{Vk, 4 1} } = [Vk, 4 + 1 ] = °- ( 23 ) 

It follows that, if a state \ fi) satisfies c’o4?/ ; ) = solV*)) then we have cr^^felV’) — — 'Sojfcl'0)- 
Analogously, if 4+llV’) — Sil'C’), we have = SiT]k\fi). Therefore, if \ip) belongs 

to sector (s 0 ,s.l+i), then r)k\if) belongs to sector (—so,sl+i)- 

To conclude, we should determine to which sectors the two degenerate ground states 
belong. If J 0 and Jl are both positive, we expect the ground-state configurations to be 
ordered. Therefore, if on one boundary the spin is directed in the positive x direction, 
we expect the same to occur at the other end of the chain. Therefore, we conclude 
that the ground states belong to sectors (+1,+1) and (-1,-1). This identification is 
supported by exact diagonalization. 

We can now classify the states. We only consider the states with s 0 = +1, to avoid 
the double degeneracy. Then, if J 0 , Jl > 0 we have: 

1) In the sector s 0 = 1, sl+i = 1, the lowest energy state is the ground state of the 
Hamiltonian and all states are obtained as ry ki .. .ry k 10), fcj > 1, with m even. The 
first excited state is rjlr}l\0) and the energy gap is A = Si + £ 2 . 

2) I 11 the sector s 0 = 1, sl+i = — 1, the lowest energy state is the first excited state 

77 J |0) of the Hamiltonian H e and all states are obtained as ry k ... 10), k t > 1, 

with m odd. In particular, the first excited state in the sector is 10), so that the 
energy gap is A = £ 2 — £i- 

The first case is relevant when considering Ci,Cl > 0 in the boundary Hamiltonian H \,, 
which are identified with Jo and Jl, respectively. The second case is relevant when 4 l 
and (l have opposite signs: the corresponding gap is obtained by identifying 4 = J 0 
and Cl = - Jl- 


3.2. Exact residts for the energy gap 

The energy spectrum of Hamiltonian H e can be computed exactly, as discussed in 


Appendix A Here we will focus on the case Jo = Jl = C, which allows us to derive the 




Quantum Ising chains with boundary fields 


7 


spectrum for parallel or opposite equal magnetic fields at the boundary, cf. Eqs. (J3J) and 
(HJ), respectively. Neglecting the irrelevant zero eigenvalue, we obtain L + 1 elementary 
excitations whose energy £ m is given by 


£ m = 2\/l + g 2 — 2g cos k ri 


(24) 


(25) 


where k m are the L + 1 solutions of the equation 

[(1 + g 2 )(1 — 2(fi) + C 4 — 2<7 cos k(l — C 2 ) 2 ] sin kL = 
g{ 1 — 2 g cos k + g 2 — C 4 ) sin [k(L + 1)]. 

We should consider real solutions in [0,7r[ (correspondingly |1 — g\ < £ m < 1 + g), purely 
imaginary solutions k = ih with h > 0 (correspondingly £ in < |1 — g\), and solutions 
k — 7T + ih (correspondingly £ m > 1 + g). Equation (I25]i has also a spurious solution 
for k = 0 for all values of the parameters, which should be discarded. For g = 1 — ( 2 , 
there is also a second solution with k = 0, which corresponds to a true excitation of 
the model. For ( = 0 we obtain the equation appropriate for open boundary conditions 

(OBC) 0, 

smk(L + l} = 1 (26) 

sin kL g 

and 1 + g 2 — 2gcosk = 0, which implies an additional zero mode for H e . This is not 
surprising as the spectrum has a fourfold degeneracy for Jo = Jl = 0: the four sectors 
are equivalent in the absence of boundary fields. 

The real solutions of Eq. (1251) are obtained by solving 


g sin k(l + g 2 — C 4 2g cos k ) 


(1 — C 2 ) 2 + g 2 ( 1 — 2( 2 ) — g[g 2 — 1 + (2 — £ 2 ) 2 ] cos k + 2 g 2 cos 2 k 


tan kL = 


The imaginary solutions are obtained by setting k = ih. Eq. (125|) then becomes 


(27) 


e 2hL = 


e h (e h g — l)[e fe (l — C 2 ) ~ 9? 
(g-e h )(e h g- 1 + C 2 ) 2 


(28) 


The corresponding energies £ = 2 + g 2 — 2g cosh h are always smaller than those 

corresponding to real momenta, hence the localized states are the most important ones 
in the determination of the low-energy spectrum. 


4. The quantum paramagnetic case 

For g > 1 (paramagnetic phase) the lowest energy state is a localized state, i.e., a 
solution of Eq. (I28p with h > 0. For large values of L we must consider the poles in the 
right-hand side of the equation. They are 

1-C 2 


e h = g, 


e h = 


9 


(29) 
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For g > 1, the second equation does not have any solution. The relevant pole is the first 
one, so that the relevant solution can be written as e h = g + e. For large values of L, 
substituting this expression in Eq. fl28j), we obtain 

CW-1) 


e = 


(g 2 - i + C 2 ) 2 


g 


Correspondingly, we find 


£ _ O CV-P q -L 

& i ^ t g 


(30) 


(31) 


<9 2 -l + C 2 )' 

The energy £\ vanishes with exponentially small corrections. The energy of the second 
excited state corresponds to the lowest value k min of k satisfying Eq. (1271) . We obtain 
fcmin = vr /L + 0(L~ 2 ) and 


£2 — 2 (g — 1) + 


gTT 2 1 


+ 0(L~ 3 ). 


(32) 


g-lL 2 

The system is therefore gapped, and the same gap is obtained for both PBF and OBF 
modulo exponentially small corrections, since £ 2 — £\ ~ £2 + £\ + 0(g~ L )\ 


A = 2(s-l) + 


gir 2 1 


+ 0{L~ 6 ) 


(33) 


g-iL 2 

Note that the same result would have been obtained for OBC, confirming that boundary 
conditions and/or boundary fields are irrelevant in the paramagnetic quantum phase. 


5. Critical-point behavior 

5.1. Finite-size scaling for g = 1 

At the critical point g — 1, there are no localized solutions, hence the two lowest energies 
£ n are obtained by considering Eq. (1271) . The equation becomes 

cot (/c/2) (</ 4 — 4 sin 2 k/2) 


tan kL = 


4C 2 — C 4 ~ 4 sin 2 k/2 


(34) 


The lowest-energy solutions correspond to momenta that scale as k ~ 1/L. Therefore, 
we can expand the right-hand side in powers of k, obtaining 


tanfci = -(hT + 0(*;), 


(35) 


which shows that tan kL must diverge for L —> 00 . This fixes k n fa (n — 1/2)7 t/L with 
n — 1,2,... Including the corrections, we obtain for the two lowest momenta 


so that 


J n tt c 2 - 4 

1 2L L 2 4C 2 ’ 

c TT 7T (fi 2 — 4 
1 = L + L2 2(2 ’ 


37t 3tt ( 2 — 4 
2 = 2L + lF 4C 2 ’ 

c 37t 37t ( 2 — 4 

2 = T + L2 2C 2 ' 


(36) 


( 37 ) 
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These results imply that the PBF and OBF gaps are 


4tt 

ApBF = ~T 


. 2tt 

Aqbf — —r 


2tt C 2 - 4 
7T — 4 


(38) 

(39) 


respectively. Note that the amplitude of the 1/L term differs from that obtained in the 
OBC case (corresponding to £ = 0) [5], 


Aqbc — — + 0(L 2 ), 


(40) 


indicating that £ = 0 is a surface critical point. 

As we have discussed in Ref. p[IJ in the OBC case, the corrections of order 1/L to 
the leading behavior can be interpreted as the effect of a nonlinear scaling field associated 
with L: the leading irrelevant boundary operator gives corrections that scale as L~ 2 . 
The same holds in the presence of boundary fields. Indeed, we can define a new length 
scale L e ff so that 

Ajr 2tT 

Apbf = j -f 0(L~z), A 0 bf = j -f (41) 

n e ff n e fr 


Using the previously reported expressions we obtain 

1 C 2 - 4 

Lett = L - - ^ 2 • (42) 

Note that the same rescaling applies both to PBF and OBF, a result which is not obvious 
in the original formulation, but which has a natural explanation in the present setting 
in which we consider H e . Indeed, the rescaling should cancel the leading correction in 
the whole low-energy spectrum of H e , hence in all low-energy excitation energies £k- 
Therefore, it should apply to both PBF and OBF. 

Since ( = 0 is a surface critical point and the boundary term is a relevant 
perturbation, we can study the crossover behavior close to ( = 0. The relevant scaling 
variable is = ( L 1 ^ 2 . Therefore, in the limit £ —y 0, L —> oo at fixed ^ we have 


A(C, L) = Lf b (( b ), (43) 

where the function fb{(b) depends on the type, OBF or PBF, of boundary fields. In the 
scaling limit defined above, momenta scale as k = h/L. Using Eq. (1341) . we obtain an 
implicit equation for b: 

i y 4 _ 7,2 

Therefore, to compute the scaling function, one determines the two lowest values of b 
satisfying Eq. (144)1 and uses them to determine the gap as before. The scaling curves 
are reported in Fig. |4] 
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Figure 1 . (Color online) Crossover curves for LA/ it as a function of C,b = l/L 1 ^ 2 . For 
(b ->0we recover the OBC case and LA/n — »• 1 . For 0 — > +oo we obtain LA/tt 2 
for OBF and LA/w — > 4 for PBF. 


5.2. Scaling behavior close to the critical point 

Let us now consider the scaling behavior in a neighborhood of the critical point g — 1. 
The corresponding scaling variable is 

w = (g- 1 )L. (45) 

The relevant equations are obtained from Eqs. ({27j) and (1281) by expanding the right- 
hand sides in the limit g —> 1, L — > oo, k —> 0 at fixed w and kL. As expected, as long 
as £ 7 ^ 0, we obtain a result that is independent of (: the FSS behavior only depends on 
the nature of the boundary conditions, but not on the specific values of the boundary 
fields. 

A localized state exists only for w > 1. If we set h — S/L, the parameter 6 is a 
solution of the equation 

5 — w tanh 5. (46) 

The corresponding excitation energy is 

Sloe = jVw' 2 - SS (47) 

Note that 6 —> 0 for w —> 1, so that S\ oc L —>■ 2 in the limit. On the other hand, for 
w —> oo we obtain 5 -A w(l — 2e~ 2w ) so that £\ oc L ~ Awe~ w , consistently with Eq. (T3Tjl . 
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W 

Figure 2. (Color online) FSS curves for LA/tt, for PBF, OBF, and OBC, as a function 
of w. For w = 0, LA/tt = 4, 2,1 in the three cases, respectively. 


The propagating modes can be written as k = e n /L, where e n are the positive 
solutions of the equation 

e n — wtane n , (48) 

with (n — 2>/2)n < e n < (n — l/ 2 ) 7 r (the solution with n — 1 only exists for 0 < w < 1). 
The corresponding excitation energies are then S n = 2 yhn 2 + e 2 . Using these results we 
obtain for the energy gap in the FSS limit: 

LA(w) = 2\jw 2 + 63 ± 2\jw 2 + e 2 for w < 0, 

LA(w) = 2 \Jw 2 + e 2 2 ± 2 \Jw 2 + e\ for 0 < w < 1 , 

LA{w) = 2 \Jw 2 + e| ± 2\/w 2 — S 2 for w > 1. (49) 

The two signs refer to the case of parallel (plus sign) and opposite (minus sign) boundary 
holds. The curves are reported in Fig. [2] together with the corresponding OBC scaling 
function. In the paramagnetic phase all curve have a similar shape — they coincide 
asymptotically — while in the phase g < 1 (w < 0), the behavior is clearly different. For 
w —> — 00 , we have LApbf(w) ~ 4|w|, LAobf(w) ~ 37 t 2 /|w|, and LAobc(^) ~ 4\w\e w . 

It is interesting to discuss the scaling corrections, that can be easily shown to decay 
as 1/L, as it occurs for ( = 0. In Ref. [2T| we argued that such corrections are not 
associated with irrelevant boundary operators, which give rise to corrections of order 
L~ 2 , but that they should be interpreted as due to nonlinear scaling fields. The same 
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holds for C 7^ 0- First, one should consider the effective length L e g defined in Eq. 

Then, one should consider the nonlinear scaling field u^ associated with g — 1 and, 
finally, rescale the energies with the sound velocity c(g). The nonlinear scaling fields do 
not depend on the boundary conditions, hence we can use the expressions reported in 
Ref. [21]: 

c{g) = 2^g u ll = (g-l)/ y /g. (50) 

Then, we define the new scaling variable 

w = Uf.LeS (51) 

and express e n and 5 as functions of w. A somewhat lengthy calculation shows that 
c(g)LA/2 can be written as in Eq. (l49jh with w replacing w, with corrections that scale 
as L~q. As expected, the leading corrections cancel out. 


6. The magnetized phase 


Let us now consider the case g < 1. In the infinite-volume limit, there exist two 
degenerate ground states |+) and |—), which differ by the value of the magnetization [5] 

m± = hm /l _ K) ±hm L _ i . 0 o(c r l 1) ) = ±m 0 , m 0 = (1 - g 2 ) 1/8 , (52) 

where h is a bulk magnetic field applied along the x direction. For a chain of finite size 
L, the degeneracy is lifted. In the OBC case, the energy difference between the two 
lowest-energy states vanishes exponentially as L increases [5j: 

Aqbc = E 1 -E 0 = 2(1 - g 2 )g L [l + 0(g 2L )] . (53) 

Moreover, m(x) =0 by symmetry, and 


G(L/2,x) 

ml 


->• 1 


for L —> oo. 


(54) 


These are the standard features of a quantum ferromagnetic phase with a spontaneously 
broken Z 2 symmetry. 

We now extend the discussion to the case in which boundary fields are present, 
considering equal PBF and OBF. To identify localized states, we consider Eq. ([28|h The 
relevant solutions are associated with the poles in the right-hand side of the equation, 
see Eq. f[29]l . It is immediate to verify that there are no poles for (fi > 1 — g, while a pole 
exists in the opposite case. Therefore, we distinguish three different cases, depending 
on whether ( is smaller, larger, or equal to 


Cc(s) = (i - s) 1/2 . 


(55) 
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6.1.1. Low-field regime, ( < (fig). For ( < (fig), the two lowest-energy excitations 
corresponding to £\ and £2 are localized. We have k = ih ( h > 0) with 

1 _ r- 

e h = s(l + 5), s =-, (56) 

9 

where s > 1 and 5 is a correction term. Let us first obtain a large-L expression for 6 
that is uniform as ( —> 0. This requires particular care as the right-hand side of Eq. (|28|) 
has a different singular behavior for £ = 0 and £ fi 0. For £ = 0, it has a simple pole as 
e h —> s — 1/g. On the other hand, for £ fi 0, it has a double pole for e h —> s. 

To compute 6 we set e h = s in all terms of Eq. (I28p that are not relevant for the 
singular behavior as L —> 00 , i.e., we write 


2 l = 0--ge h ) 

C l9 (e h -sfi' 

with 

r = s(s-g) 

1 g(s 2 ~ l) 2 ' 

Equation (T57T) is a second-order algebraic equation of e h , whose solutions are 

_2 l 

e h = s + (±\/ACig( 2 s 2L + g 2 - g^j . 

The energies £ 2 /4 — 1 + g 2 — 2g cosh h are correspondingly 


(57) 


(58) 


(59) 


1 , 2 = C 2 (l ~ g 2 ~ C 2 ) _ 3(1 - s 2 ) 3 s- 2 L 
4 ± (1-C 2 ) 2 (g-s)s 3 


(±^/iCxgCs 21 + g 2 - g') (60) 


and £\ = £ + , £2 = £-■ 

There are now two interesting cases. First, we take the limit L —» 00 at fixed £ fi 0. 
We can then simplify Eq. (16(1 . obtaining 



C 2 (1 


(1-C 2 ) 


a ±c ^(i 


-,-L 


g)s 3 


-Vctg- 


(61) 


These results imply 

Apbf = £2 + £\ = 4£ ^2 ^ ^ ( g L )’ (62) 

Aobf = £2 — £\ = -j- 1 - s L + 0(s 2L ). (63) 

(s - g)s 2 

In the PBF case the gap is finite and proportional to £ for small fields. Instead, Aobf 
vanishes exponentially, as in the OBC case (which corresponds to vanishing boundary 
fields). 
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A second possibility consists in considering the FSS limit around ( = 0. As 
discussed in Ref. [22], the ratio A PB f(C)/^obc is expected to become a function of 
n — £/Aobc (to compare with the formulae of Ref. [22] note that the energy associated 
with the perturbation is here proportional to £, without additional factors of L ). In 
this limit s L ( ~ g L ( is proportional to hi. Using the expression for Aobc reported in 
Eq. (l53jh we obtain 

S± 


= ^1 + 16(1- g 2 )n 2 TL 

^OBC ^ ^ 


so that 


Apbf(C) 


= Cl + 16(1 - j 2 )k 2 . 


(64) 


(65) 


^OBC 

Note that the scaling function agrees with that predicted in Ref. [22.]. We can also 
consider the ratio Aqbf(C)/^obc dr the same limit, but in this case we obtain that the 
ratio is simply equal to 1, i.e., independent of k. The reason can be easily understood. 
For OBF the energy associated with the perturbation is zero—the contributions of the 
two boundaries cancel—hence n = 0 identically. 


6.1.2. Large-field phase, Q > ( c (g)- For ( > ( c (g) there are no localized states. The 
secular equation (1271) can be rewritten as 


tan kL = s'mkf(k,(,g), /(0, C,flO 7 ^ 0. 

For L —> oo, the solutions can be parametrized as 


It follows that 


^n 

tan — « 

Lj 


/lit/ U/ n 

k = ~l + U‘- 


sin (fi +j^)f(k„, C,g). 


For large values of L we obtain finally 

a n = an = vm/(0, C, g) = nn- 


g(i-g + C 2 ) 


( 66 ) 

(67) 

( 68 ) 

(69) 


(l-g)(l-g-( 2 )' 

Note that the correction term a n diverges as ( —> ( c (g) and has a hnite limit as ( —* oo. 
Moreover, the explicit form of the Hamiltonian is only relevant for the corrections. It 
follows 

2 ag Tin 2 




1 - g L 2 1 - g L 3 


+ 0(L ~ 4 ) 


We obtain finally 


Apbf — 4(1 — g) + ——— zr^ + 0(L 3 ), 


Aqbf — 


1 — g L 2 

3 g tt 2 + 6g 2 (l - g + C 2 )vr 2 + 


(70) 


(71) 

(72) 


1 -g L 2 (l — 5 , ) 2 (i ~ g — C 2 )E 3 
For OBF, the gap vanishes for L —> oo. However, while the approach is exponential for 
C < Cc(fl'), we have A OB f L 2 for large fields. For PBF, the system is gapped as for 
( < £ c . However, in this case corrections are of order L~ 2 . 
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6.1.3. Intermediate state, ( = (fig). The gap equation is particularly simple. We 
obtain a solution cos k — 1 and 


sin(L + 1 )k 
sin Lk g 

Hence £\ = 2(1 — g) with no L dependence. As for the other states we obtain 

gm 2 7T 2 1 2gm 2 'K 2 1 


(73) 


£m +1 — 2(1 — g) + 


1-gV (1 -g) 2 L^ 


(74) 


with m = 1,... L. It follows A PB f = 4(1 — g) + 0(L~ 2 ), as for ( > ( c . However, 
corrections, of order L~ 2 , are a factor of five smaller compared to the large-field case. 
For OBF we obtain the asymptotic behavior 


A - 9 ^ 

‘Aqbf — --— 


2 gn 2 


+ 0(L~fi. 


(75) 


1 - 9 L 2 (1 - g) 2 L 3 

As in the large-field case, Aobf scales as L~ 2 . However, the prefactor differs by a factor 
of 3 from that obtained for ( > ( c , cf. Eq. (1721) . 


6.2. Magnetization and correlation function 

We wish now to compute the ground-state magnetization and the corresponding two- 
point spin-spin correlation function in the magnetized phase, i.e., for g < 1. In the 
PBF case, the system is magnetized, hence we expect 777(7) and G(i,j) to be both 
independent of i and j, except possibly close to the boundaries, and to be equal to m 0 
and /Up, respectively, where m 0 is given in Eq. (I52j) . Thus, the boundary magnetic fields 
have the only role of lifting the degeneracy. 

The OBF case is more interesting. In Appendix B we have fully characterized the 
ground state of the model for small values of g. We find that for it is ferromagnetic for 
( < ( c (this is consistent with the results for the gap discussed above). For ( > ( c kink 
states [I] are the relevant low-energy excitations, while for ( = ( c the ground state is 
a superposition of kink and ferromagnetic states. The exact knowledge of the ground 
state allows us to compute the ground-state magnetization and the correlation function 
perturbatively in g. For OBF we find that 777(7) = (cy* ) vanishes in the low-field case 
except at the boundaries, depends linearly on i for ( = (fig), and varies in a simple 
fashion for ( > (fig). We expect the same space dependence for all finite values of 
g < 1 in the large-L universal limit, modulo a multiplicative renormalization. If we set 
x = i — L/2 and I = L/ 2, so that the center of the chain corresponds to x — 0, we 
predict 


m{x) 

777 0 

mfx) 

777 0 

m(x) 

7?7o 


x 1 irx 

= - + - Sill — 

I 7T t 

X 


= 0 


C > Cc, 
C = Cc, 
C<Cc 


(76) 
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Figure 3. Plots of m(x)/m o versus x/t for OBF in the ferromagnetic phase g < 1. 
Here L = 21 + 1, x = 0 at the center of the chain. Top: Results for g = 1/2, 
C = 1 and g = 3/4, £ = 3/4 (large-field kink phase). Middle: Results for g = 3/4 at 
C = Cc = \/l — g ■ Bottom: Results for g = 3/4 and 5 = 1/4 (small-field magnetized 
phase). The dashed lines barely visible on top of the data correspond to the theoretical 
curves GSD. 
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These relation should hold in the limit x —* oo, £ —* oo at fixed x/£, except possibly at 
the boundary, i.e., for xji= ±1. Here m 0 is the bulk magnetization of the system in the 
infinite-volume limit defined in Eq. (T52T) . It is interesting to note that the linear behavior 
of the magnetization for ( = ( c is due to the fact that the ground state is translationally 
invariant: the space probability of the kinks is independent of x. Due to the presence 
of the boundary fields, the translation invariance of the kink space distribution is not 
trivial and, indeed, it does not hold for ( > ( c . 

To check these predictions and study the approach to the asymptotic behavior, we 
perform numerical simulations using the density matrix renormalization-group (DMRG) 
method For convenience, we consider chains with odd L, setting L = 2£ + 1 (since 
we are considering ferromagnetic interactions, the low-energy properties do not depend 
on this feature, unlike the antiferromagnetic case, see, e.g., Ref. M). The results for 
the magnetization are shown in Fig. [3l For ( < ( c , the magnetization vanishes in a 
large interval around x = 0, assuming positive and negative values only close to the 
boundaries. Note that the region in which m(x) 0 shrinks as L increases, conforming 
the validity of Eq. (176|) in the large-L limit at fixed x/£. For ( > ( c we show results 
for two values of g. Once we renormalize m(x) by using mo, i.e., we consider the ratio 
m(r)/m 0 , results fall one on top of the other, confirming universality. The results are 
in full agreement with the expressions (l76lh 

The same arguments apply to the correlation function G(i, j). Extending the small- 
g results of|Appcndix B|to the whole low-g phase, we predict 


G(i,j) , 1, M ,-m 

-o— = 1 - \m(i) — m(j)\ 


mn 


m 0 


(77) 


in the large-L limit. Therefore, the correlation function should be constant in the 
magnetized phase, G(i,j ) = rn ^, it should depend linearly on \i — j\ for ( = ( c , while it 
should have a sinusoidal dependence in the kink phase. Again, we check these predictions 
by comparing them with DMRG data. They are reported in Fig. [4j The DMRG 
data approach the theoretical predictions with increasing L, confirming the theoretical 
predictions. 


7. Magnet-to-kink transition for 9 < 1 

As already discussed in Ref. [23] . in the case of OBF and for g < 1, an interesting 
universal scaling behavior is observed for ( —> ( c (g)- Such a transition is uniquely 
characterized by the nature of two phases occurring for C < Cc and C > Cc : on one side 
the ground state is ferromagnetic, while on the other side kink states are the relevant 
low-energy excitations. This characterization is confirmed by the analysis of the Ising 
chain on a ring with a bond defect: at the boundary of the two phases one observes [23] 
a behavior analogous to the one discussed here. Here, we will report the computation of 
the energy gap— the result was already presented in Ref. [23]—and additional results 
for the local magnetization, correlation function, and entanglement entropy. Moreover, 
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x/1 


Figure 4. Plots of G{y = 0,x)/mo versus x/£ for OBF (L = 21 + 1, x = y = 0 at 
the center of the domain). Top: Results for g = 3/4, £ = 3/4 (large-field kink phase). 
Middle: Results for g = 3/4 at ( = ( c = y/1 — g. Bottom: Results for g = 3/4 and 
C = 1/4 (small-field magnetized phase). The dashed lines barely visible on top of the 
data correspond to the theoretical curves obtained by using Eqs. OS) and dzzD. 
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we will discuss the density of the fermionic modes in the equivalent fermionic picture 
of the model |28] [29], showing that the transition can be interpreted as a localization 
transition of the fermionic states at the boundaries. 

The crossover behavior around the transition is parametrized in terms of the scaling 


variable [23] 


C, = ^^£(C-Co). 

(78) 

Then, in the large-L limit we have 


Aqbf(C) ~ A 0 bf((c)/a(Cs), i/L PH Cs)« 

(79) 


In general, the scaling functions /a(Cs) and ffifis) are expected to be universal modulo 
a (/-dependent normalization (which is independent of the observable considered) of 
the argument ( s . The g dependent prefactor appearing in Eq. d78|) is the required 
nonuniversal factor. It is determined in Sec. 17.11 analyzing the scaling behavior of 
the gap Aobf- According to renormalization-group theory, this normalization factor 
is independent of the observable considered, hence also the scaling functions of other 
observables should be independent of g apart from a multiplicative factor (not present if 
renormalization-group invariant ratios are considered), once ( s is used as scaling variable. 
The results we will present for the magnetization and the two-point correlation function 
confirm this prediction. 

Beside the gap Aobf(C)> we also consider the integrated correlation x defined in 
Eq. (jHJ). It satisfies the scaling relation 

X = rn 2 0 Lf x (Q, (80) 

where f x (Cs) is universal, i.e., g independent. These results allow us to define the RG 
dimensions of the different operators. In general, we expect 

A 0 Bp(C) = i-‘/i[i' I (C~Cc)], X = A x L d -y’J x [LX((-Q], (81) 

Comparison with our results allows us to identify z — 2, — 1, and y a = 0. The 

latter result is consistent with the fact that a bulk magnetic held coupled to cA 1 - 1 drives 
a first-order quantum transition, see, e.g., Ref. [22]. 

By defining ( s as in Eq. (1781) we have implicitly assumed that we are fixing g and 
varying ( close to the critical value Cc(d)- In some contexts, however, it is more natural 
to consider a different point of view. The boundary held strength ( is hxed (it should 
satisfy ( < 1) and g is varied in a neighborhood of the critical value g c = 1 — ( 2 . In this 
case we have ^ 

C = 2(i-c 2 ) ^ ~ 9 ^ L ' ( 82 ) 

where the ( dependent prefactor guarantees that scaling functions are the same for all 
values of (. 
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Figure 5. (Color online) Scaling behavior of the gap Aobf around C = £ c . The 
DMRG estimates (points) of Aobf(C)/^OBf(Cc) approach the universal theoretical 
curve for L —> oo (full line). 


7.1. Scaling function for the energy gap 

To determine the scaling behavior of the energy gap Aobf(C)j we first simplify the 
secular equation ([25]) . by taking the limit C —» Cc at fixed Cs- In this limit k scales as 
1 /L, so that we set k — z/L. For L —> oo the secular equation becomes 

4 QgZ + (4^s — z 2 ) ta ,nz = 0. (83) 

Expressing tan z as a function of tanz/2, we hnd that the solutions of Eq. (l83j) satisfy 
one of the two equations 

tan | = -2^’ (84) 

tan - — -A (85) 

If z a i < z a 2 < ... and Zbi < Zb 2 < Zb 3 ■ ■ ■ are the positive solutions of Eqs. (184|) and (f85l) . 
respectively, it is easy to verify that Zbi < z al < Zb 2 < ■ ■ ■ for ( s > 0 and ( s < —1. For 
— 1 < ( s < 0 we have instead z a i < Zbi < z a 2 ■■ ■■ For ( s < 0 there are also localized 
solutions with k = iu/L. The parameter u satisfies one of the two equations 

tanh- = - —, (86) 

tanh— = ——. (87) 

2 u 

Eq. (I87|) has a positive solution Ub for ( s < 0, while Eq. (1861) has a positive solution 
u a satisfying u a < Ub only for ( s < —1. These results allow us to determine the 
scaling function defined in Eq. (1791) . We obtain /a(C) = (-al - Z bl)/^ 2 for Cs > 0, 
/a(Cs) = (zl 1 + u 2 b )/TT 2 for -1 < Cs < 0 and f A (( s ) = (■ u 2 - u 2 a )/TT 2 for Cs < -1- 
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Figure 6. (Color online) Rescaled local magnetization m(x)/mo versus x/t for 
C s = -10,-3,0,3,10. 


We can easily determine the asymptotic behaviors. For small ( s we have 

/a(Cs) — Id y + 0(Gs)> (88) 

7T Z 

while for ( s —> oo we have 

/a(C) = 3-^ + 0(C 2 )- (89) 

Finally, for ( s —> — oo we have 

/a(C) « ^C 2 e 2Cs - (90) 

7T z 

The theoretical curve is shown in Fig. [5] The DMRG data confirm the theoretical 
calculation: the numerical estimates of the ratio Aobf(C)/^obf(Cc) clearly approach 
the scaling curve /a(Cs) with increasing L. 


1.2. Magnetization and two-point function 


Let us now consider the behavior of the correlations close to ( c . The local magnetization 
m(i) and the correlation function G(i,j) are computed exactly in the limit g —> 0 in 


Appendix B.4 As we have already done in Sec. 16.21 we extend these expressions to all 
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o 

c. 


Figure 7. (Color online) Crossover scaling curve for x/(toqL) versus ( s . We report 
the theoretical prediction (dashed line) and DMRG data for g = 1/2 and g = 3/4 
to confirm universality. Note the limiting values: f x ( 0) = 1/2, f x (— oo) = 1, and 
f x ( oo) ~ 0.297, indicated by the horizontal dashed lines. 



c. 


Figure 8. (Color online) Crossover scaling curve £/L = fi(Cs) for the correlation 
length £. We report the theoretical prediction for the asymptotic scaling behavior and 
DMRG data for g = 1/2 and g = 3/4 to confirm universality. Note the limiting values: 
fd 0) = l/y/4S, ffi- oo) = 1/^24, and / € (oo) w 0.0985. 
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values of g satisfying g < 1, by replacing m{i) with m(i)/m 0 and G(i,j ) with G(i,j)/mfi 
Therefore, we predict for ( s > 0 

m(x) z 

mo z + sin z 

where z is the smallest positive solution {z = Zbi) of Eq. (j85j) , x = i — L/2, and £ = L/2, 
as before. For ( s < 0 we have instead 

m{x) u 

mo u + sinh u 

where u is the solution {u = ufi) of Eq. (j57|) . As for the correlation function G(i,j ) 
it satisfies Eq. f|77|) . ft is important to stress that the previous expressions, although 
derived for small values of g, are expected to be exact for any g < 1. For instance, since 
m(x)/mo is a dimensionless renormalization-group invariant ratio, it should scale as 

m(xX,g,L) = m 0 / m [a(g)« - QL,x/L], (93) 


1 , 

rux\ 

X 

— smh 


T ~7 

u 

\ £ J 

£ 


(92) 


1 . 

(zx\ 

, x 

- sm 

— 

~\~ — 

z 

^ £ ) 

£ 


(91) 


where f m is a universal, hence g independent, function. The specific features of the 
model enter in this expression only through the normalization nonuniversal factor a(g) 
and the critical value (fig). The factor a(g) does not depend on the observable, hence 
it can be determined by using any quantity. We use the gap Aobf, for which we are 
able to obtain exact results for any g < 1, defining ( s = a(g)(( — (fi)L. Therefore, f m as 
a function of ( s and x/L is g independent. Therefore, our results, obtained in the limit 
g —> 0, hold for any value of g satisfying g < 1. 

The previous expressions for rri(x) hold in the limit x —* oo, L —* oo at fixed xji 
and smoothly interpolate between the expressions valid in the different phases. Note 
that we have m{%) —> ±mo for xjl —>■ ±1 in all cases. This result does not however 
necessarily apply to the boundary points, i.e., for x/£ = ±1, as the magnetization is 
not continuous at the boundary. For instance, in the magnetized phase, m; = 0 in all 
interior points, while mi and rrii = —mi are nonuniversal functions of the field strength 

c 

A graph of m{x)/mo for several values of ( s as a function of xjl is reported in 
Fig. [6j As expected, for ( s —» — oo the magnetization approaches zero exponentially in 
an interval that is centered in the middle of the chain and that widens as |£ s | increases. 
More precisely, for ( s —» — oo we have u ~ — 2( s and 


m(x) = 


sinh 2( s x/£ 


(94) 


sinh 2 ( s 

which shows that rn(x) vanishes except in two tiny intervals close to the boundaries of 
width £/(21 C s |) - 

To verify the prediction (177]) for G(i,j ) we compute the scaling functions associated 
with x and (. They are reported in Appendix B.4 and compared with DMRG data in 
Figs. [7] and El Results with different values of g fall on top of the theoretical curve when 
plotted versus ( s . This confirms our predictions for G(i,j ) and the magnetization, as 
well as the correctness of the nonuniversal prefactor appearing in the definition of ( s . 
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Figure 9. (Color online) Scaling of the subtracted fermionic density at the center of 
the chain. We plot DMRG data of Ln s { 0) for g = 3/4 and g = 1/2 versus ( s . The 
full line corresponds the asymptotic analytic curve, cf. Eqs. (1981) and (l99l) . Note that 
Ln s ( 0) approaches 2 for ( s —>■ oo, 0 for Q —> — oo, while it is equal to 1 for ( s = 0. 


7.3. Equivalent fermionic picture 

It is interesting to reinterpret our results in the equivalent fermionic picture of the model. 
In the ferromagnetic phase, i.e., for ( < ( c , the lowest eigenstates are superpositions of 
Majorana fermionic states localized at the boundaries p8j ,29]. In finite systems, their 
overlap does not vanish, giving rise to the splitting A ~ e ~ L / l o _ The coherence length l 0 
diverges at the kink-to-magnet transitions as ~ | In s| ~ £ c — £, a behavior analogous 
to that observed at the order-disorder transition g — > 1~ where Iq 1 I l n S' I • 

Within this picture, it is natural to consider 


n{i) 


i + (d 3) > 


(95) 


which can interpreted as the particle density in the fermionic representation of the Ising 
chain, see Eq. m- Its large-L limit is given by [5] 


1 1 r .. g + cos k 

nh = —I-/ dk — , 

2 2 tt J o + g 2 + 2g cos k 

We consider the difference 

n s (x) = n(i) - n h , 


(96) 

(97) 


where x — i — L/ 2, as before. Such a quantity has been computed for small values of g 
in Appendix B.4| The universality argument we have used for m{i ) and G(i,j) can be 
applied to n s (x). We therefore predict 


zx\ 

1 + COS-] 


Ln s {x) 


z + sin z 


(98) 
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Figure 10. (Color online) Subtracted fermionic density Ln s (x ) versus x/t for 
( s = —10, —3,0, 3,10 (a; = 0 at the center of the chain). 

for ( s > 0 and 

Ln s (x ) =--— (l + cosh —(99) 

u + sinh u V t) K J 

for Q s < 0, where z and u are determined as in the case of the local magnetization. 

To verify this prediction we compute n s ( 0) numerically for g — 1/2 and g = 3/4. The 

numerical results are compared with theory in Fig. 0 The agreement confirms the 

analytic prediction. 

The subtracted fermionic density Ln s {x ) is reported in Fig. [TO] as a function of 
x/t. For positive values of ( s it has a maximum at the origin and shows that fermions 
are delocalized. For ( s = 0, Ln s { 0) is constant, while for negative values of ( s , there 
is a significant enhancement at the boundary. For ( s —> —oo, as in the case of the 
magnetization, Ln s (x) is different from 0 only close to the boundary, in an interval of 
width ^/2|£ s |. It is interesting to consider the behavior of Ln s (x ) at the boundaries, 
i.e., for x —» ±£, which is shown in Fig. [TT] In the kink phase, such a quantity scales 
as 7 t 2 /(2(//) for ( s — » oo. On the other hand it diverges as — 2£ s for ( s — » — oo, that is 
when the magnetized phase is approached. As such, it represents a physically relevant 
order parameter for the transition, which distinguishes the two different phases. 
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- 4-2 0 2 4 6 8 10 

Cs 

Figure 11. (Color online) Subtracted fermionic density Lnb = Ln s (±£) at the 
boundaries of the chain versus ( s . We have Lnb (( s ) ~ 7r 2 /(2^) for ( s —► oo, 
Lnb(Cs) ~ —2Cs for ( s —> —oo, and Lrib(O) = 1. 



Figure 12. (Color online) Scaling of the subtracted half-chain entanglement entropy. 
The DMRG data appear to approach an asymptotic scaling curve. The scaling 
curves for different values of g match after a global nonuniversal rescaling: fixing 
c(g = 3/4) = 1 the optimal matching is obtained for c(g = 1/2) ss 1.85. 
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Finally, we discuss the behavior of the half-chain entanglement entropy [32] - We divide 
the chain into two connected parts [—^, 0] and [1,£] and consider the von Neumann 
entropy 

S = -Tr P[-e, 0 ] In P[_y 0 ], (100) 

where 0 ] is the reduced density matrix of the subsystem [—£, 0], i.e., 

P[-t, 0] = Tl '[iy] p, (101) 


and p is the density matrix of the ground state. The entanglement entropy S approaches 
a constant in the large-L limit, except for the critical case g = 1 where it increases 
logarithmically [32]. Its large-L limit for OBC is given by flHJ 


S'oBc(fi') — Y2 


In 


160 ^ 


+ 1 


I{z) = 


dx 


In Fig. 


0 \/(l — x 2 )(l — Z 2 X 2 ) 
we plot the difference 


g 2 \ 4J(g)J(yT^7) 

2 ) TT 


( 102 ) 


S, = S((,g,L)-S 0BC (g), 


(103) 


which shows the scaling behavior 


Ss^L-'fsiCs), (104) 

where fs is a universal function apart from a multiplicative nonuniversal constant. 
These results imply that the large-L limit of the half-chain entanglement entropy does 
not depend on the boundary conditions, being the same for OBC and OBF. However, 
the O^L^ 1 ) corrections are affected by the boundary term (l2lh giving rise to a nontrivial 
universal behavior around £ c . 


7.5. Comparison with the existing results for two-dimensional classical models 

The results we have obtained allow us to derive exact scaling functions for the classical 
two-dimensional Ising model in a strip Lx 00 in the presence of two equal magnetic 
fields on the boundaries. This model has been extensively studied, especially in the 
case of opposite magnetic fields. In this case, the wetting (or interface localization- 
delocalization) transition, which is the analog of the magnet-to-kink transition we are 
considering here, has been extensively studied [331 EH E3 [36], [371 EH ESI US] EH E2] • We 
can therefore check some of our results, comparing them with the existing ones in the 
literature. 

The results for the gap Aobf obtained in Sec. 17.11 can be used to derive the scaling 
behavior of the longitudinal correlation length defined by = 1/A 0 bf- Such a 
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quantity was computed in the restrictive solid-on-solid model in Ref. [37]. Their result 
is completely consistent with ours. If our scaling quantities z and u defined in Sec. 17.11 
are identified with the quantities y/2Go, appearing in Ref. we obtain exactly 

the same expression for £||/L 2 . 

Several papers study the magnetization profile, which gives information on the 
behavior of the interface between the two coexisting phases. The behavior in the large- 
held phase was determined in Ref. [35] by using the restrictive solid-on-solid model. 
The result is in agreement with our prediction ( 1731 ) for ( > £ c . The linear behavior 
observed at ( = ( c was also derived at the wetting transition by using the solid-on- 
solid model [39]. No exact results are instead available for the magnetization profiles 
in the whole crossover region, which, in the classical case, is usually parametrized 
in terms of [T — T w (h)\L, where h is the boundary held and T w the corresponding 
wetting temperature (in the quantum-to-classical mapping g corresponds to T, hence 
this scaling corresponds to varying g at hxed £). However, the numerical results for 
the magnetization reported in Refs. [40], jTl, 42] (and, in particular, Fig. 2 of Ref. [40]) 
appear to be in full agreement with our prediction shown in Fig. [6] 

8. Conclusions 

We consider the one-dimensional quantum Ising chain in the presence of a transverse 
magnetic held g [5] and of boundary magnetic helds aligned with the order-parameter 
spin operator. In particular, we consider helds with equal strength (, which have either 
the same (PBF) or the opposite (OBF) direction, cf. Eqs. (J3]) and (SJ). We assume the 
chain to have length L and derive the fini te-size behavior of several quantities in the 
limit L —» oo as a function of (. 

We derive analytic predictions for the gap in all phases as a function of the boundary 
held strength. In the paramagnetic quantum phase, the leading behavior is independent 
of the boundary conditions. At the critical transition g — 1 , finite-size scaling depends 
on the boundary fields. The point ( = 0, corresponding to open boundary conditions, 
is a surface critical point. We derive analytic expressions for the scaling functions that 
parametrize the crossover behavior for g — y 1 and any ( (bulk critical behavior), and 
for g = 1 and ( > 0 (surface critical behavior). In the quantum ferromagnetic phase 

g < 1, if the boundary fields are oppositely aligned, the finite-size behavior drastically 
changes as £ varies [23]. For small £, the system is in a ferromagnetic phase in which the 
gap decreases exponentially with the system size. On the other hand, for large fields, 
kink propagating modes of momenta 1/L are the relevant low-energy excitations, so 
that the gap decreases as 1/L 2 . The two phases are separated by a critical transition 
at ( = £ c (g ). Close to this transition, low-energy properties show a universal scaling 
behavior in terms of the scaling variable L(£ — £ c ). The transition is only characterized 
by the nature of the coexisting phases. Indeed, the same transition occurs in Ising rings 
in the presence of a localized bond defect [23] . 

It is interesting to interpret our results in the context of two-dimensional classical 
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Ising models, using the quantum-to-classical mapping. As is well known pQ, the quantum 
Ising chain of length L corresponds to a classical Ising model on a strip of width 
L , the transverse field playing the role of the temperature. Therefore, our results 
provide exact predictions for the Ising model in an L x oo geometry with boundary 
fields. Our results for the gap at the critical point can be used to predict the FSS 
behavior of the longitudinal exponential correlation function = 1/A as a function of 
w = A(T — T C )L and of Q — B^L 1 ^ 2 , where A and B are nonuniversal model-dependent 
constants. Moreover, our extensive results at the magnet-to-kink transition for g < 1 
for OBF provide predictions for the finite-size scaling functions at the classical wetting 
transition (sometimes also called interface localization/delocalization transition [2?]) in 
two dimensions. In particular, we have exact predictions for the FSS behavior of the 
magnetization profile and of the correlation function of the layer magnetization in the 
whole crossover region close to the transition. 

In the context of the wetting transition, one may also consider a second important 
surface parameter, which is usually referred to as surface enhancement |33, .21] • While a 
boundary magnetic field is the most relevant magnetic perturbation that breaks the 
Z 2 symmetry at the boundary, surface enhancement is the most relevant (energy¬ 
like) perturbation that is Z 2 invariant. In the context of the quantum chain, such 
an interaction can be mimicked by adding a term (/qcr® + h^o ®), localized at the 
boundary, to the general Ising-chain Hamiltonian. It is not difficult to generalize the 
results presented here to this more general case (results for the spectrum in the absence 
of boundary magnetic fields appear in Refs. mm)- At the critical point, this would 
allow us to study the behavior at the ordinary and extraordinary surface transitions (in 
the absence of surface enhancement the point g — 1 , £ = 0 corresponds to the so-called 
special transition 031). Moreover, a complete understanding of the phase behavior for 
g < 1 in the absence of a bulk magnetic held should be possible. In the quantum case, 
these results would also be relevant for Ising quantum rings with localized site defects, 
a distinct case with respect to that considered in Ref. [ 23] . 


Appendix A. Spectrum determination 


Appendix A.l. Preliminary matrix residts 

In this appendix we collect results on the spectrum of some matrices that are relevant 
for the discussion. We begin by considering the n x n matrix (we write it for n — 5) 

(d b 0 0 0\ 


A n (d, b ) = 


d 

b 


b 
d 

0 0 b 

\0 0 0 


0 0 
b 0 


d 

b 


b 

d J 


(A.l) 


To determine its eigenvalues, we compute 

a n (d, b) = det (A n (d, b) — XI). 


(A.2) 
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a n = 2b cos fca n _i — fe 2 a n _ 2 . 


(A.3) 


The solution is 

a n b~ n = Cl e ikn + c 2 e~ ikn , (A.4) 

where C\ and c 2 are arbitrary constant. If we require cq = 2b cos k and a 2 = 
5 2 (4cos 2 k — 1), we obtain 


xe 


ik 


Cl = 


It follows 


(In b 


2 sin k 

sin k(n + 1) 


c 2 = c v 


sin k 


(A.5) 


(A.6) 


The eigenvalues are solutions of the equation a n = 0. We obtain n eigenvalues given by 

nm 


A m = d — 2b cos k r 


k 


n + l’ 


(A.7) 


with m — 1,n. 

Let us now compute the eigenvectors. Let (aq,...,a n ) be the eigenvector 
corresponding to eigenvalue A m . It satisfies the relations 


2ai cos k m + a 2 = 0, 

aj-i + 2aj cos k m + a 3 +i =0, j — 2,... n — 1 
c^ 7 i—i T 2a n cos k m 0. 


The solution of the recursion relation (1 A.9 11 is 


(A.8) 

(A.9) 

(A.10) 


aj = (-l) j ( ei e ik ™ j + e 2 e~ ikmj ), 


(A.ll) 


where e\ and e 2 are two arbitrary constants. Requiring a 2 
e\ + e 2 = 0, so that 

otj = (—l)' 7 e 1 sin k m j, 


—2cosfc m ai, we obtain 


(A.12) 


where we have rescaled e\ for convenience. Requiring the eigenvector to have unit norm, 
we obtain 

OLj = (-1) J \ — 2 — sin k m j. (A.13) 

V n + 1 

We also determine the eigenvalues of the matrix 


A n (d, b, e) 


(e b 0 0 0\ 

b d b 0 0 

0 b d b 0 

0 0 b d b 

^0 0 0 b ej 


(A.14) 
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which differs from A only for two matrix elements: 
Setting again 

a n = det (A n - XI), 


A n n and A U)nn are equal to e. 


(A.15) 


and A = d — 2b cos k, we obtain 

a n = (e — d + 2bcosk) 2 a n _2 ~ 26 2 (e — d + 26cos k)a n _ 3 + 6 4 a n _ 4 . (A. 16) 

Using the recursion relation (1A.3I) . we obtain finally 

a n — (1 — 5 2 )a n + 2b5(6 cos k — l)a n _i, (A. 17) 


where 5 = (d — e)/b. 


Appendix A.2. Secular equation 

To determine the spectrum of Hamiltonian (191) . we compute the eigenvalues of the matrix 
C defined in Eq. (1201) . We write them as 

A = g 2 + 1 — 2g cos k. (A. 18) 


Note that nothing forbids A to be larger than |1 + g\ or smaller that |1 — g\, hence cos k 
is not necessarily bounded between —1 and 1. Hence, we must consider real values of k 
with 0 < k < n and complex values of the form k = ih and k = it + ih with h > 0. 

We are not interested in the zero eigenvalue of C, hence we will only consider the 
(L + 1) x (L + 1) matrix formed by the first (L + 1) rows and columns, which we label 
with the symbol Cl+i- Then, we introduce a matrix Dl+i defined as (here L — 4) 


Dl+i — 


(g 2 +1 g 0 

g g 2 +1 g 

0 g g 2 + 1 

o o g 

\ 0 0 0 


0 

0 

g 

g 2 +1 

g 


The corresponding secular equations are 

c n = det {C n — XI) d n = det (D n — XI). 


0 

0 

0 

g 


\ 


g 2 + ij 


The determinant d n has already been computed in Appendix A.l 


d n = g‘ 


t sin k(n + 1) 
sin k 


(A.19) 


(A.20) 


(A.21) 


To determine c n , we expand the determinant with respect to the first and the last rows. 
We obtain the following recursion relation: 

c n = (g 2 + 1 — Jo ~ 2 g cos k)( 1 — J 2 L — 2g cos k)d n - 2 

+ g 2 (g 2 + 1 - JqJI - 2g cos k - 2gj$ cos fc)d„_ 3 + g 4 J 0 2 dn- 4 - 


(A.22) 
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A simpler expression is obtain by eliminating d n _ 4 and d n _ 3 , and expressing the result 
in terms of d n _\ and 2 . This can be done by using the recursion relation satished by 
the determinants d n , see Eq. (1A.3h with b = g. We obtain finally 


c n = - (1 + g - J 0 J L - 2gcosk)d n _ 1 


(A.23) 


+ [(1 + g 2 )( 1 - J 0 2 - Jl) + JlJ\ - 2g(l - J|)( 1 - Jl) cos k] d n _ 2 . 

Note that the recursion relation (1A.23D is symmetric with respect to the exchange of J 0 
and Jl- The eigenvalues of C n are obtained by requiring c n = 0. Setting J 0 = Jl = C and 
n = L + 1, we obtain Eq. f[25lh Note that we have simplihed sin A; in the denominators, 
hence Eq. (l25l) has a spurious solution for k = 0. A true solution with k — 0 occurs 
only when the coefficients of d n _\ and d n _ 2 in Eq. (1A.23|) both vanish for k — 0. This 
only occurs when J% — = 1 — g. 


Appendix B. Perturbative analysis in the small-g limit for opposite 
boundary conditions 

The low-energy behavior of the model with Hamiltonian © can be understood 
analytically in the limit g —> 0. For g — 0, it is trivial to compute the spectrum of 
Hamiltonian (OQ) in the presence of the magnetic boundary term (121) , since H is diagonal 
in the basis in which is diagonal. Assuming Ci = — Cl = — C an d C > 0, there are 
two family of states that control the low-energy spectrum. First, we should consider the 
magnetized states 

I+> = |1,1,1,•••,!), 

|-) = |-i,-i,-i, (b.i) 

then, we should also consider the (L — 1) states (we call them kink states p]) 

|i>* =l-i, i,i, ■■■,1), 

|2)fc = I - l.-l.l.1), 

\L — l)fc = | — 1, —1,...,—1,1). (B.2) 

If Hq is the Hamiltonian for g = 0 we have (J = 1) 

H 0 \±) = —(L - 1)|±>, = -{L- 3 + 2C)|i) fc - (B.3) 

The value of ( determines which of the states is the ground state of the system. For 
( < 1, the ground state is doubly degenerate and spanned by |+) and |—), while for 
c > 1 the (L — 1) kink states are the lowest-energy ones. For £ — 1, the magnetized 
and the kink states are degenerate. We wish now to determine how this picture changes 
when the perturbation 

H„ = sY. a ? } < B4 ) 


is added. 
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Let us first understand the behavior for ( < 1. For g — 0 the ground state is doubly 
degenerate. Such a degeneracy is lifted by perturbation (1B.4j) . Note, however, that 
(±|P™|+) vanishes for any n < L. Therefore, the energy gap is proportional to g L , in 
agreement with the exact results obtained in Sec. | 6 l 

Let us now consider the corrections to the ground state which are proportional to 
g. At order g, |+) mixes with the states generated by H g |+), so that we can write 

■0+ — a \+) + fi\ — 1,1,1,.. ■, l)i + <S|1,1, ■ ■ •, 1, —1 )l + 

L—l 

+ 7 ^|l,l,...,-l,...,l,l) i , (B.5) 

i=2 

where 7 , 5 are of order g and state |.. .)* is dehned so that afi .. .)* = |.. .)* for 
i j and < 7 j| .. .)* = — | .. .)*. Requiring fi> + to be normalized, we immediately obtain 
a = 1 + 0(g 2 ). The coefficients 7 , and 5 are fixed by the eigenvalue condition 

H„\+) = (Eo - ff„)W+ - |+)). (B.6) 


We obtain 


P = 

The analysis for |- 


--- 7 = - 5 =--- (B 7) 

2(1-C)’ 4’ 2(1 + 0' ' 

-} is analogous and leads to = TP z fi> + . Once the degeneracy is 


lifted, the ground state should be an eigenstate of TP,, hence it should be of the form 


1 

7 ! 


(V+ ±^-)- 


(B. 8 ) 


We are now in the position to compute correlation functions on the ground state. For 
the local magnetization, we obtain 

mi = 0 for i = 2,... L — 1 

»i\ = —m L = S 2 - /3 2 = - 9 _ ^ . (B.9) 

This result shows the the ground state does not show a local magnetization, except 
at the boundaries. For i 1, T, the local magnetization m[i) vanishes as in the zero- 

held case. The applied held ( is too small to destabilize the low-energy state. We can 

analogously compute the correlation function G(i,j): 

G(l,L) = 1 ,2/3= - 2<5 2 = 1 — (B.10) 

G(i, L) = G(l, t) = 1 - f - « 2 - 2 7 2 = 1 - ( 5 (B.11) 
G(i,j) = 1 -4y 2 = 1 - g 2 /A = m 2 0 + 0(g 3 ), (B.12) 

where 1 < i j < L. Except on the boundaries, the correlation function is equal to the 
square of the bulk magnetization mg, as in the zero-held case. 
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Let us now consider the large-field region in which the ground state is (L — 1) degenerate 
and the relevant states are the kink states (IB.2D . The degeneracy is lifted when H g is 
included. The full Hamiltonian restricted to the subspace spanned by the kink states 
has the form (we write it for L — 5) 


H = 


^0 

~9 

0 

0 \ 

-.9 

E 0 

-g 

0 

0 

~9 

Eq 

-9 

1° 

0 

~9 

Eq ) 


(B.13) 


where E 0 = — (L — 3) — 2(. The eigenvalues and eigenvectors of this matrix have been 

(B.14) 


computed in |Appendix A.l| There are L — 1 energy states 

E m — Eq 


7T m 

2q cos ——, 
L 


with m = 1,... L — 1. The ground state corresponds to m — L — 1. The corresponding 
eigenfunctions are 


3 =1 


-l) 2 sin k m j\j)k, 


(B.15) 


with k m = nm/L. Using this expression, we can compute the average of cr,- 1 ^. For the 
average on state m, we obtain 


(m\a^\m) = — (m\a'i > \m) = — 1, 


Wi™\ - 


(B.16) 


( m \ a i \ m ) = sin2 k mj - y- Sin 2 k m j , 

j<i j>i 

where 2 < j < L — 1. If we now consider the ground state, in the large-L limit, we can 
rewrite Eq. (1B.16D as 

2j 1 . 2vr j 

rrij = — 1 + —-sm ——. 

L it L 

If we define x = i — L/2 and l — Lj 2 we obtain the more symmetric form 

. . X 1 . 7TX 

m (^) = 7 + - sm —■ 

Analogously, we can compute the correlation function. We obtain 
{m\a\ a L \m) = —1, 


(B.17) 

(B.18) 

(B.19) 


/ I (1) (1)I \ / I l 1 ! I \ 

{rn\a\ crj \m) = \m), 


(!) I 


/ | ( 1 ) ( 1 ) | \ / | | \ 
{m\a L a j \ m) = (m|crj | m), 


(!) I 


(m 


o-j'Vy \m) = j sin 2 k m i - j ^ sin2 ^ + j ^ sin 2 fc m z, 


^<,7 j<i<k i>k 

where 2<j<k<L — 2. Focusing again on the ground state and taking the limit 
L —> oo, if x = j — L/2, y = k — L/2 , we obtain 


G(j, fc) = 1 


\x-y\ 

1 

7 TX 

sin —— 

Try 

— sin — 

t 

7T 

£ 

£ 


= 1 — |m(j) — m(/c)|. (B.20) 
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We should now discuss the behavior of the system in the intermediate case in which 
kink states and ferromagnetic states are degenerate. For g — 0 this case corresponds to 
C = 1. As discussed in Sec. |6l when the g dependent term is added, the intermediate 
case corresponds to ( c = Q1 — g ~ 1 — g/2. Therefore, to discuss the intermediate case, 
we set ( = 1 — g/2 and consider the model restricted to the (L + l) subspace spanned 
by 1°) = |-), I i)k (k = 1, • • •, L - 1), | L) = |+). If E 0 = — (L - 3) - 2( c is the energy of 
the kink states, the Hamiltonian can be written as (here L — 4) 


( Eq - g 

-9 

0 

0 

0 

\ 

~9 

E 0 

-9 

0 

0 


0 

~9 

Eo 

-9 

0 


0 

0 

-9 

Eq 

-9 


V o 

0 

0 

-9 

Eo- 

V 


(B.21) 


The eigenvalues can be computed using the results of |Appcndix A.l| Indeed, H has the 
same form as matrix A defined in Eq. flA.141) . Since 5 = —g, the secular equation (1A.17I) 
becomes 

sin k(L + 1) 


(1 + cos k)- 


sin k 


= 0, 


which implies 


k = 


Tim 


m — 1,... L + 1. 


(B.22) 


(B.23) 


L + l 

As the energy of each mode is given by E = E 0 + 2g cos k, the ground state is obtained 
by taking k — ix. The ground-state energy is Eq — 2g and the corresponding eigenvector 
is simply 

r 

l 


I GS) = 


+ 1 j=0 




(B.24) 


It is easy to compute the local magnetization and the correlation function. We find 

2 i 


m(i ) = — 1 + 


(B.25) 


L + l ’ 

which shows that m(i ) varies linearly between 1 — 2/(L + 1) and —1 + 2/{L + 1). If we 
set, as usual, x = i — L/2 and t = L/2 we obtain in the large-L limit 


x 


m(i) = -. 

The two-point correlation function is also easily computed 


(B.26) 


G(i,j) 


2| i-j 
L + l 


1 — |m(i) — m(j) |. 


(B.27) 
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Let us now study the crossover behavior. We consider again the same basis as in 

\ 


Appendix B.3 The Hamiltonian becomes (here L = 4) 

(E 0 


H = 


9 + 2(?C s/L 


-9 


\ 


0 

0 

0 


(B.28) 


0 
0 
0 

-9 

Eq — g + ZgCs/L/ 

The secular equation is obtained by using the results of |Appcndix A.l[ Eq. (1A.17I) . We 
set b = —g, d = E 0 , e = E 0 — g + 2 g( s /L, and E = E 0 + 2g cos k. Since 5 = 2 ( S /L — 1, 
we obtain 


-g 

Eq 

~9 

0 

0 


0 

~9 

Eq 

-9 
0 


0 

0 

-9 

Eq 

9 


f" 


c 


2C, 


-j- ) sin A; (L+2)+2 ^1 — -j- 


1+ 1 


2C 1 i- 
-j- I cos k 


sinA;(L+l) = 0. (B.29) 


In the kink phase in which ( s > 0, k varies between k = n for ( s = 0 and k = tt — 7r/ L 
for ( s — > oo. Hence, we write k = tt — z/L. Expanding the secular equation to order 
1/L 2 we obtain 

4 ( s z + (4<C S 2 — z 2 ) tan z = 0, (B.30) 

which coincides with Eq. (j83]h The solutions of this equation either solve tan z/2 = 2 ( s /z 
or tan z/2 = —z/{2(, s ). The ground state corresponds to the lowest value of z that is a 
solution of Eq. (IB.301) . hence it satisfies tan z/2 = 2( s /z. Note the limiting values: 

2 ~ 2a/G for Cs ->■ 0, 

z ~ 7r(l — 1/C) for C ~^ +oo. (B.31) 


We wish now to compute the ground-state eigenfunction, which we express as 
(«i,..., qll+i) in the basis | j), j = 0,... L, defined in Appendix B.3 If k = n — k = z/L, 
the coefficients satisfy the relations 

ai(2cosfc — 1 + 2( S /L) — a 2 = 0, (B.32) 

ctj-i — 2aj cos k + a.j +1 = 0, (B.33) 

ctL — oi L+i (2 cos k — IT 2</ s / E) = 0. (B.34) 

The solution of Eq. (IB. 3 31) is 


aj = e\ cos kj + e 2 sin kj. 


(B.35) 


The constants e± and e 2 are hxed by condition (IB.321) and by the normalization condition 
£«? = !' In the large-L limit these two conditions give 


ei = Za 


(4C 2 + 4 Cs + z*)L' 


e 2 = 2( s 


(4C/ + 4 Cs + zQL' 


(B.36) 

(B.37) 
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We can now compute m[i ) and G(i,j). For the local magnetization we have 


m(i) = a\- a? L+1 + a ) ~ 


a. 


1 <j<i i<j<L 

In the large-L limit, using tanz/2 = 2 (, s /z, we obtain 

... 0 


m[i) = 


1 zx x 
~z sm T + l 
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(B.38) 


(B.39) 


z + sm z 

where x = i — L/ 2, i = L/ 2. Note that the prefactor guarantees that m(i) is equal to 
±1 at the two boundaries. We can also compute the correlation function, which satisfies 
the relation 

G(i,j) = 1 - | m(i) - m(j)\. 

Correspondingly, we obtain 

X z 2 — 2 + 2 cos z + 2z sin z 
L 

and 


e 

L 2 


2 z(z + sin z) 
z 4 + 24 + 12 (z 2 — 2) cos z + 4 z(z 2 — 6) sin 0 


(B.40) 


(B.41) 


(B.42) 


A8z 2 (z 2 — 2 + 2 cos z + 2z sin z) 

Let us finally consider n s {i) defined in Eq. ([97j). For g —>■ 0, nh — 1/2, so that 

/o\ 

n s (i) = {a\ ) /2. Using the definition we obtain 

n s (i) = OLiOti+x. (B.43) 

Substituting the explicit expression of cq and using tanz/2 = 2 ( s /z, we find 


Ln s (i) = 


1 + cos 


zx\ 

TJ 


(B.44) 


(z + sin z) 

The above-reported results apply to the kink phase ( s > 0. For ( s < 0 the ground state 
is a localized state with k = it + iu/L, where u satisfies the equation 

4 ( s u + (4(/ 2 + u 2 ) tanh u = 0. (B.45) 

The solutions of this equation satisfy either tanhu/2 = —2 ( s /u or tanh u/2 = —u/(2( s ). 
Since E = E 0 — 2g coshu/ L, the ground state is obtained by considering the largest 
positive solution of the equation. It is then easy to prove that it satisfies tanh u/2 = 
—2 (s/u. The solution for ( s < 0 is then obtained by analytic continuation of the results 
obtained for ( s > 0. It is enough to replace z with —iu. We therefore obtain 

u 


m{i) = 
Ln s (i) = 


x 1 ux 

. ~ n + - smh — 

u + smh u \ t u t 


u 


1 + cosh 


Note that u 


u + sinh u 
2( s for ( s —>• —oo, so that 

sinh2^ s x/£ 


ux\ 

T ) 


m(i) = 


sinh 2( s 


Ln s (i) = 


2( s cosh 2( s x/£ 
sinh 2( s 


(B.46) 

(B.47) 

(B.48) 


in this limit. 
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